Forward-backward multiplicity correlations in relativistic heavy-ion collisions in a 

superposition approach 
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We analyze the multiplicity correlations between distant forward and backward rapidity regions in 
relativistic heavy-ion collisions in a superposition framework, where the particle production occurs 
through independent emission from correlated sources. This allows for inferring information on the 
long-range forward-backward correlations of the sources in the earliest phase of the collision, based 
on the experimental information on the statistical features of the observed particle distributions. Our 
three-stage study incorporates the fluctuations of parton production in the early phase, the effect 
of intermediate hydrodynamic evolution, as well as fluctuations in the production of particles at 
freeze-out. We investigate the dependence of the results on the features of the overlaid distributions 
and hydrodynamics, as well as the centrality dependence. We confirm the previous conclusions of 
the substantial forward-backward event-by-event fluctuation of sources implied by the data from the 
STAR Collaboration. Predictions based of the Glauber model of the initial phase for the forward- 
backward multiplicity correlations in Pb+Pb collisions at the LHC are made. 

PACS numbers: 25.75.-q, 25.75Gz, 25.75.Ld 
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I. INTRODUCTION 

The forward-backward (FB) multiplicity correlations 
in relativistic heavy-ion collisions, recently measured at 
RHIC [J-13]) were followed with a number of theoreti- 
cal studies [5rfl9j. The primary goal of these analy- 
ses is to get insight into the space-time dynamics of 
the earlies stages of the reaction, probed via the long- 
range correlations. A simplified statistical understand- 
ing [1, 0, [l4[ of the problem is gained when one considers 
initially formed sources (wounded nucleons [20l-|22j| . pos- 
sibly amended with binary collisions [23, 24], Glasma |25l - 
[28j . dual strings [H[29j]), which later lead to particle pro- 
duction occurring independently (and with the same dis- 
tribution) from each source. Such superposition models 
lead to simple relations between the statistical features 
(moments, correlations) of the distributions of sources 
and of the produced particles, which involve properties 
of the overlaid distribution as parameters [3] ■ 

In this paper we analyze in detail the predictions of 
this framework, tailoring it closely to the popular de- 
scription of the relativistic heavy-ion dynamics based on 
three stages: initial phase, hydrodynamics, and statisti- 
cal hadronization (for a review see, e.g., 30]). The key 
assumption is that the emission from a source is uni- 
versal, i.e., occurs from the same statistical distribution 
independently of the centrality class of the collision. This 
leads to very simple formulas for the statistical measures, 
in particular for the correlation coefficient between the 
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numbers of particles produced in the forward and back- 
ward rapidity bins. Parameters in these relations depend 
on the features of the overlaid statistical distributions 
and the properties of hydrodynamics. 

We use the derived relations in two ways. First, using 
the experimental data from the STAR Collaboration 0- 
4] for the correlation and the scaled variance, we obtain 
an expression for the correlation of sources in the early 
phase, involving one free parameter dependent on the un- 
known features of the overlaid distributions and hydrody- 
namics. This straightforwardly generalizes the treatment 
of Brogueira and Dias de Deus [5[ and Bialas, Bzdak, and 
Zalewski [lij , who use a Poisson distribution in the two- 
step approach. We confirm the conclusions of Ref. [l4| 
of the possible large hidden event-by-event asymmetry in 
the distribution of sources in the RHIC setup. We carry 
out our analysis for the Au+Au and Cu+Cu collisions at 
,,/snn = 200 GeV, where the data are available. 

Second, we use the relations in the opposite direction, 
starting from the properties of the source distribution 
and computing the FB correlation of the produced par- 
ticles. We use the Glauber framework for this purpose. 
Predictions, involving one free parameter, for the Pb+Pb 
collisions at ^/snn — 2.76 TeV at the Large Hadron Col- 
lider (LHC) are made that way. We predict a decrease 
of the FB correlations of particles with centrality, which 
follows from the decrease of the scaled variance of the 
sources. The future data from the LHC will verify this 
scenario. 

The outline of the paper is as follows: In section (jn} we 
briefly recall the three-stage approach, consisting of the 
early production, the hydrodynamic evolution, and the 
final freezeout. We then derive the basic relations link- 
ing the statistical features of the sources and of the final 



2 



particle distributions. In Sec. pil[) we use the data from 
the STAR Collaboration for the Au+Au and Cu+Cu col- 
lisions to obtain the FB multiplicity correlations of the 
sources. Finally, in Sec. (|IV[) we present our predictions 
for the FB correlations of particle multiplicities to be 
experimentally analyzed with the LHC for the Pb+Pb 
system. We use the approach introduced in Sec. (JTTJ) 
and the distributions of sources obtained with the mixed 
model [HSl from GLISSANDO [p. 



II. THREE-STAGE APPROACH 

Much of success in the description of the relativistic 
heavy-ion dynamics has been achieved in a three-stage 
model (for a review see, e.g., consisting of 



1. 



2. 



Early production, modeled in terms of the Glauber 
approa ch [ 20l -|24l. l32| . Glasma [25-28], string forma- 
tion [jl.l29l|. etc. 

Hydrodynamic evolution, starting from the initial 
conditi on p rovided by stage 1 (for reviews see, e.g., 
30l 3S. I3J . and for the recent event- by- event stud- 



ies 



mm) 



3. 



Statistical hadronization, carried out at freezeout 
right after the hydrodynamic phase ends, e.g., [451 - 



The approach leads to a proper description of such 
measured quantities as multiplicities, spectra, harmonic 
flow coefficients, femtoscopic properties, etc., J57i], thus 
is viewed as a practical framework. Statistical fluctua- 
tions are generated in phases 1 and 3, while phase 2 is 
assumed to be deterministic (although one may include 
fluctuations also in this phase [58j|). Of course, the most 
interesting are the fluctuations in the initial stage, as 
they refer to the important physics at the earliest times 
and may help to discriminate between various physics 
approaches, while the fluctuations at hadronization form 
a "statistical noise" which should be disposed of. While 
each of the stages is physically involved, including nu- 
merous physical parameters, and needless to say, takes 
a huge effort to simulate numerically, certain statistical 
aspects, as we shall see, can be understood and classified 
in rather simple terms, displaying the possible scenarios 
in the FB correlations of the earliest phase. 

Throughout this paper we use the generic concept of 
sources, which may be viewed statistically as the density 
of partons (or fields) in the initial phase, turning into 
the entropy density of the fluid cell in the hydrodynamic 
phase, which in the end, at freeze-out, gives birth to the 
observed hadrons streaming to the detectors. In the con- 
text of the FB multiplicity fluctuations, an important 
assumption is the sufficient kinematic separation of the 
forward (F) and backward (B) rapidity windows, such 
that the particles produced in from the F source do not 
fall into the B window, and vice versa. This allows us 



to trace the evolution of a cell with fluctuating sources. 
Suppose in the early-production phase we separate the F 
and B cells in spatial rapidity, with the original number 
of sources s denoted as sf and sb, respectively. These 
numbers fluctuate event-by-event and may be correlated, 
which is what we eventually want to assess. In the initial 
production mechanism these sources produce partons p, 
whose density supplies the initial condition for hydrody- 
namics. 

Assuming that the production occurs from each source 
in the same manner, we then have 



PA 



A = F.B, 



(1) 



where the random variable fj,i is the number of partons 
produced from the source i. As mentioned, we assume 
that the distribution of //$ is universal, i.e., does not de- 
pend on the location of the cell, and that the produc- 
tion from different cells is independent from one another. 
Then the formulas for the superposition model follow (see 
Appendix |A"|) : 



(Pa) = (v)(sa), 
vai(p A ) = v&r(fj,)(s A ) + (p) 2 va.r(s A ), 
cov(pf,Pb) = (/i) 2 cov(s F , s B ), 



(2) 



where A = F, B. 

The sources p yield the entropy density which consti- 
tutes an event-by-event initial condition for the collec- 
tive evolution via hydrodynamics. The hydrodynamic 
evolution is dynamically complicated, however, it is de- 
terministic. Thus the evolution of the cell with initially 
p sources yields h sources at freeze-out (p and h may be 
thought of as entropy contained in the cell), where h is 
a function of p0 The phase-space location of the fluid 
cells is evolved due to the hydrodynamic push, however, 
the separation between the far-lying F and B regions is 
maintained. If the event-by-event fluctuations are not 
too large, we may expand 



to(p)+h( P -(p)) + 0((p-(p))- 2 



(3) 



where f, are parameters depending on dynamical prop- 
erties of hydrodynamics. The higher-order terms may 
be dropped if p is sufficiently close to (p) . The constant 
term is written in the form to (p) . Then 



(h A ) = toipA), 
vax(h A ) = £?var(pA), 
cov(h F , h B ) = t\cov(p F ,p B )- 



(4) 



1 Strictly speaking, the hydrodynamic evolution depends not only 
on the number of sources p, but also on their spatial distribu- 
tion, hence the same values of p may lead to somewhat different 
h. Such nuances can only be included in a fully numerical simu- 
lation. We do not expect them to be relevant for our analysis. 
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Equation Q requires some discussion. If the response 
of hydrodynamics is such that h ~ t, i.e. the entropy 
in the considered cell after the evolution is proportional 
to the initial entropy, then to ~ t\. The success of hy- 
drodynamics in reproducing the rapidity spectra starting 
from an initial Glauber condition (see, e.g., (59|) supports 
this scenario. Moreover, it suggests that to and t\ are 
independent of centrality in the range of applicability 
of hydrodynamics. Secondly, since in viscous hydrody- 
namics the entropy is produced, we have in a given cell 
(h) > (p), which implies t Q > 1. The production of en- 
tropy depends on the properties of hydrodynamics, but 
is not very large [5!| (25-50%), such that we expect to 
to be only somewhat larger than 1. To summarize, we 
arrive at the estimates 



to — t\ 
to > 1. 



1, 



(•5) 



The parameters to and t\ should be treated as specific 
to a given reaction or rapidity bin, but independent of 
centrality. 

Finally, statistical hadronization is carried out at 
freezeout. We assume that a given cell emits n hadrons 
into a region of phase-space with some statistical distri- 
bution superimposed over h. Similarly to the analogous 
mechanism in the initial phase, each of the h sources 
emits independently m hadrons with the same distribu- 
tion. In reality some mixing may occur and particles may 
be emitted from different cells into the same kinematic 
region. However, if the F and B cells are well separated, 
this effect is negligible. Thus we have 



n A = y^m i; 



(6) 



and 



(n A ) = (m)(h A ), 
v&r(n A ) = v&r(m)(h A } + (rn) 2 va,v(h A ). (7) 
cov(nF,n-s) = (m) 2 cov(hF,hB)- 

The three-stage model and the notation introduced above 
may be summarized with the following diagram: 



init. production hydro . hadronization 

s — > p — > h — > n 



Joining Eqs. (HO yields 

(n A ) = a(s A ), 
var(ru) = /3(s A ) +7var(s j4 ), 
cov(np, Ub) = 7Cov(sj?, sb)- 

where the combinations of constants arc 



(8) 



It is convenient to introduce the scaled variance 
uj(x a ) — vax(x A )/{x A ) and the correlation coefficients 
p(xf,xb) — cov(xf,xb)/{(j(xf)o'(xb))- Then, for sym- 
metrically separated bins where <j(xf) = &(xb), we can 
write the relations 



u(n A ) 
p(n F ,n B ) 



5 + kui(sa), 
p(sf,s b ) 
1 + X/uj(s a )'- 



(10) 



with 



S = /3/ct = oj(m) H (m)w(/Lt), 

to 

K = y/ot = —(n){m), 
to 



(11) 



A = /?/7 



tooj(m) w(/x) 



Formulas (|8I X0[) express the statistical properties of the 
event-by-event distributions of the produced particles in 
the F and B bins via the properties of the distribution of 
the original sources. Note that the relation between the 
correlation coefficients depends on a single combination 
of the unknown parameters of the overlaid distributions 
and hydrodynamics represented in A, which makes the 
qualitative analysis simple^ 

We may also write the reversed formulas expressing the 
properties of the initial sources through the properties of 
the final particle distributions: 



(sa) = ~(n A ), 

a 

w(sa) = ^ 

K 

p(n F ,n B ) 



(12) 



p(s F ,s B ) 



1 - 5/w{n A )' 



A feature following from the assumptions in our anal- 
ysis is the independence of the parameters a, /3, 7, 8, k, 
and A on the centrality class. This dependence resides 
entirely in the statistical properties of the distributions 
of s or n. 

Finally, we note that the inclusion of finite experimen- 
tal acceptance a amounts to overlaying yet another sta- 
tistical distribution over n, with mean a and variance 
a(l — a). As a result, the parameters a, j3, and 7 are 
changed, however, the form of Eqs. (|51 [T2"|) remains unal- 
tered. 

The basic methodology is as follows: Eq. ([12"]) involves 
three independent parameters: a, k, and A (note that 
5 = kX). Thus, knowing from the experiment (n A ), 
cj(n A ), and p{np,nB) and from the model (s A ), oj(s a ), 
and p(sf,sb) (at a given centrality) allows us to solve 



a = t Q (p)(m}, 

(i = to(p)v&r(m) +t l (m) var(/x), 
y = tl(p) 2 (m) 2 . 



(9) This feature is a derivative of the form of Eq. J3}. Note also that 

overlaying more distributions preserves the form of Eq. (181 10ft . 
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FIG. 1. (Color online) Forward-backward correlation of the 
initial sources in Au+Au collisions inferred from the STAR 
data [3,01, plotted as a function of centrality for several values 
of the 5 parameter in Eq. (|12[) . The case 5 = (joined with 
the line) corresponds to the STAR data for p(rif,ng). 



the equations for a, k, and A. Comparing the values at 
various centralities serves as a consistency check0 For 
this program to be feasible, one needs the complete ex- 
perimental data involving (ua), w(tia), and p{n F ,riB) 
at all centralities, which, unfortunately, is not the case of 
Rcfs. 3, 4], where u)(tia) is provided only for two central- 
ity classes. For that reason the above program is carried 
out only partially in the following Sections. On the model 
side, we need (sa), w(sa), and p(s f ,sb)- 



III. ANALYSIS OF THE STAR DATA 

We now use Eq. (JT2J) and the data from the STAR 
Collaboration @,H| for the Au+Au and Cu+Cu collisions 
at RHIC at the energy ^/snn = 200 GeV to compute the 
correlation of sources p(s f ,sb)- We use here the data 



3 More appropriately, one should simultaneously solve the equa- 
tions at all centralities in the \ 2 sense. 
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FIG. 2. (Color online) Same as Fig. [T] for the Cu+Cu colli- 
sions. 



for the the well-separated symmetric bins at Ai] = 1 
and A77 = 1.8. As already mentioned, the experimental 
variance var(s^), necessary for the analysis in Eq. (|12l) . 
is provided only for two centrality bins: 0-10% and 40- 
50%. The values of (sa) at various centralities are read 
off from Fig. 4.1 of [|. 

The results are shown in Figs. [1] and [U where the 
FB correlations of the initial sources in the Au+Au and 
Cu+Cu collisions is plotted as a function of centrality for 
several values of the S parameter of Eq. (fP2)l . The case 
(5 = corresponds to the STAR data for p(n F ,nB)- Nat- 
urally, increasing 5 leads to larger p(s f ,sb)- We cannot 
use too large values of 5, as p(s F , Ss) < 1 and, as follows 
from Eq. (Tig), 



5 < uj{n A )[l - p{n F ,n B )}. 



(13) 



Numerically, it gives the maximum allowed values for S 
at c = — 10% equal to 1.6 and 1.3 for Au+Au and 
Cu+Cu, respectively, while at c = 40 — 50% the maxi- 
mum values are only slightly larger than 1. Insisting on 
the independence of 8 on centrality, we find (from the 
STAR data for c = 40 - 50%) that 6 ~ 1 saturates the 
bound (PI). 

We note that the value 6=1 and the corresponding 
p(s f ,sb) — 0.8 for the Au+Au case, c = — 10%, is 
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FIG. 3. (Color online) Forward-backward correlation of the 
initial sources in Au+Au collisions obtained from Eq. (I10|) 
with various values of A. 
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FIG. 4. (Color online) Same as Fig. [3] for the Cu+Cu colli- 
sions. 



consistent with the number obtained in Ref. [141 ] . where 
(in the notation of that paper) 

£_ = ((s F - s B ) 2 ) = 1- p(s F ,s B ) _ 1 (u) 
D+ \j((s F +s B ) 2 ) V l + p(s F ,s B ) ~ 3' 1 ' 

This means that the fluctuations of the sources are quite 
substantially decorrelated, i.e., p(sf,Sb) is less than 1. 

We may now use the value of S to obtain information 
on the underlying parameters of the overlaid distribu- 
tions and hydrodynamics. From Eq. (fTTj) it follows that 

i 2 

uj(m) + -^(m)uj(p) ~ 1. The weak bound deduced from 
here is that the distribution of m cannot be wider than 
Poissonian, i.e., u(m) < 1. 

Above, we have used exclusively the experimental in- 
formation to conclude the dependence of p(sf,sb) on 
centrality, which is possible only for the two bins where 
the experimental values of uj(jia) are provided. We now 
fill the missing points by using some information from 
the Glauber model as well. For that purpose we rewrite 
Eq. (fTO)) in the form 

p(s F , s B ) = p(n F , n B ){l + X/uj(s a )), (15) 

where on the right-hand side we use the experimen- 
tal data for the correlation coefficient and the mixed 



Glauber model 31] result for uj(sa)- For the RHIC en- 
ergies yj snn — 200 GeV we take the mixing parameter 
a = 14.5% and (Jnn = 42 mb. The results of our proce- 
dure are shown in Fig. [3]|4] 

The main conclusion of this section is a gradual de- 
crease of p(s F , sb) with centrality, which, as enforced by 
the data, drops to nearly zero at c = 40 — 50%. The effect 
occurs for both the Au+Au and Cu+Cu reactions. 

We note that the first measurements of the FB mul- 
tiplicity correlations at RHIC were carried out by the 
PHOBOS collaboration [l|. Since the statistical measure 
used in that work is different form the correlation coef- 
ficient, we do not use the PHOBOS data in the present 
analysis. 



IV. PREDICTIONS FOR THE LHC 

Finally, we use Eq. (|TU1) to predict the FB correla- 
tion for the Pb+Pb collisions at the collision energy of 
s/snn — 2.76 TeV at the LHC. The values of uj(sa) are 
obtained from the mixed model Glauber implemented in 
GLISSANDO. In particular, we take the mixing param- 
eter a = 15% and onn = 67 mb. Figure ([5]) shows our 
predictions, which will be verified when the LHC data for 
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FIG. 5. (Color online) Predictions for the ratio of the FB 
multiplicity correlation of measured particles to the FB corre- 
lation of sources, plotted as a function of centrality, obtained 
from Eq. (|15[) a t different values of A. The scaled variance 
u)(sa) is taken from GLISSANDO for the mixed model. The 
simulation is for the Pb+Pb collisions at the LHC energy of 
= 2.76 TeV. 



the FB multiplicity correlations are analyzed. We note a 
graduate fall-off of p(np,ng)/p(sp,SB) with centrality. 
This fall-off is due to the fact that the scaled variance 
oj(sa) decreases with growing centrality. We note that 
the fall-off in Fig. [5] is much weaker than in the STAR 
data in Figs. 11121 in particular at c = 40 — 50% the points 
in Fig. [5] assume about a half of the central values. 



V. CONCLUSIONS 

In this paper we have extended the analysis of cor- 
relations in superposition models, previously made in 
Refs. @, I [3, to the case of the three-stage approach 
consisting of the early production, hydrodynamics, and 
statistical hadronization. Simple formulas, linking the 
statistical properties of the FB correlations in the data 
and in the original sources have been derived. The ef- 
fect of hydrodynamics may be, under reasonable assump- 
tions, incorporated in terms of just two parameters inde- 
pendent of centrality. The relations between the FB mul- 
tiplicity correlations of the produced particles and the 
initial sources involve a single parameter, which collects 
the features of the overlaid distributions and hydrody- 
namics. These one-parameter formulas allow to verify 
the model of the production early phase with the exper- 
iment, under the proviso that the data are sufficiently 
complete in providing, for each centrality, not only the 
FB correlation coefficient, but also the variance of the 
number of particles in the F and B bins. 

We have confirmed the findings of Ref. [l4[ that the 
STAR data [1 [| suggest a substantial degree of decor- 
relation of the sources in the F and B bins, which means 
large event-by-event FB fluctuations in multiplicity. A 



similar feature is seen in the Cu-Cu data. 

The awaited FB correlation analysis with the LHC 
data will shed further light on the early production mech- 
anism. The statistical method presented in this paper 
is directly applicable to that case. We argue that the 
complete data consisting of the average multiplicity and 
variance in the forward and backward rapidity bins, as 
well as the forward-backward correlation coefficient, will 
allow for a verification of production models of the early 
phase. 



This research was supported by the Polish National 
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Appendix A: Superposition model 

In this Appendix we derive the relevant statistical for- 
mulas. Let the number of produced particles n be com- 
posed of independent emissions from s sources, 

s 

n = uij. (Al) 

i=l 

Here rrii is the number of particles produced by the ith 
source from some universal distribution. Then the well- 
known relations follow: 

(n) = (s)(m), (A2) 
var(n) = (s)var(m) + (m) 2 var(s). (A3) 

We give for completeness the derivation. Introduce 

Srrii — rrii — (m), with (5rrii) = 0. (A4) 

Then 

s s 

var(n) = (V^(5m, + (m) ) (Smj + (in))) — ((s)(m)) 2 
i=i j=i 

s s s 

— (V^ <Smf) + ( SmiSnij} + 2{m) (V^ 8mi) 

s s 

+ (™) 2 (EE)-( S ) 2 H 2 - (A5) 
1=1 j=l 

The second and third term in the last equality van- 
ish due to Eq. (|A4[) . Also, from the independence of 
the production from different sources, the first term is 
equal to (s)var(m). Finally, using the obvious fact that 
Et=i E 3 S =i = s 2 we obtain Eq. (lAll) . 

Next, we look at the covariance between two well- 
separated bins, which means (mjTOj) = (w) 2 , with i and 
j belonging to two different bins. We have 

si s 2 

(nin 2 ) = (y] mi } j rrij) = (m) 2 (s 1 s 2 ), (A6) 
i=i j=i 
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and For the correlation coefficient it follows that 

p{nx,m) = , , — ■ (A8) 



cov(ni,na) = (m) cov(si,S2). (A7) 



<m)w(si) V (m)w(s 2 ) 
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